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Abstract 

We prove that lattice QCD generates the axial anomaly in the continuum 
limit under very general conditions on the lattice action, which includes the 
case of Ginsparg- Wilson fermions. The ingredients going into the proof are 
gauge invariance, locality of the Dirac operator, absence of fermion doubling, 
the general form of the lattice Ward identity, and the power counting theorem 
of Reisz. The results generalize in an obvious way to SU {N) lattice gauge 
theories. 
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1. Introduction 

In the continuum formulation of QCD (or QED) it is well known that different 
gauge invariant regularization schemes all yield the same expression for the axial 
anomaly. A non-perturbative regularization is provided by the lattice. Any candi- 
date for a lattice discretization of QED or QCD should reproduce the correct axial 
anomaly in the continuum limit. For lattice QED with Wilson fermions the axial 
anomaly has been first studied by Karsten and Smit These authors showed that 
the origin of the anomaly was an irrelevant term in the lattice Ward identity. In 
it was shown, using the small-a-expansion scheme of ref. @], that in the limit of 
vanishing lattice spacing a this naively irrelevant contribution is necessarily given by 
the D ^ A limit of the dimensionally regulated continuum triangle graph. That the 
correct anomaly is in fact generated for any lattice discretization of the QED action 
satisfying a very general set of conditions, which include the case of Ginsparg- Wilson 
fermions, has been recently shown in . It is the purpose of this paper to generalize 
this result to the non-abelian case. Although for concreteness sake we shall consider 
the case of QCD, it will be evident from our discussion that our proof goes through 
for any SU{N) gauge theory. 

The paper is organized as follows. In the following section we discuss the general 
form of the lattice axial vector Ward identity, which is one of the important ingredi- 
ents which will be needed for the proof. As we shall see, its precise structure, which 
depends on the particular discretization of the lattice action, need not be known, 
but only very general properties thereof. In section 3 we then state the theorem we 
want to prove, and discuss the consequences of the various (very general) conditions 
on the action required by the theorem. The proof of the theorem is presented in 
section 4. 

2. Ward Identity 

Consider the following general form of the lattice action for QCD 

S = SG[U] + J2'^{x){Du{x,y)+m)'ilj{x) , (1) 

x,y 

where S'G[f/] is the standard Wilson gauge field action and Du{x,y) is the Dirac 
operator (a matrix in Dirac-spin and colour space). U denotes collectively the link 
variables U^{x) = exp[igaA^{x)], with a the lattice spacing, and ip, i/j are three 
component colour Dirac fields. Ylx ~ Yin where n labels the lattice sites. The 
action is assumed to be gauge invariant, and to possess the discrete symmetries of 
the continuum theory. The Dirac operator has the following formal expansion in 
the gauge potentials 

Du{x,y) = J2 ^^£^""(^'2/1^1 ■ ■ ■ ^")^mU^i) ■ ■ ■ ^ZM . (2) 

It can always be decomposed into a chirally symmetric, and a chiral symmetry 
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breaking (sb) part as follows, 



Du{x, y) = Du{x, y)sym + Du{x, y)sb , (3) 

where ^ 

Du{x,y)sym = ^[Du,l5\'y5 , (4) 

Du{x, y)sb = \ Wu, 75)75 ■ (5) 

Any candidate for a lattice action should possess the correct continuum limit. It 
therefore follows that for a — > 0, Du{x,y)sym li-iD^lA], where -D^[^] is the co- 
variant derivative, while Du{x,y)sb vanishes in the continuum limit. 

Let 55' be the variation of the action induced by the (standard) colour singlet 
axial transformation 

Sifj^x) = ze(x)75'0(x) . 

Stfj^x) = ie{x)tl){x)'^^ . (6) 

Because the measure is invariant under this transformation, it follows that 

< 5S >u= , (7) 

where < 5S >u denotes the expectation value of SS in the external background field 
{U^{n)}. As a consequence of the Poincare Lemma on the lattice the variation 
of the action can always be written in the form 

6S = iY,uix)[-d;:jl{x) + 2mj,ix)+Aix)'\ , (8) 

X 

where ^5(2;) = '0(x)75'0(x), -0(0;) 7^75?/' (x) for a — > 0, and A(x) is an 

irrelevant operator vanishing in the continuum limit, dj^ is the (dimensioned) left 
lattice derivative. Hence the Ward identity for the colour singlet axial vector current 
will be of the general form 

< dl:jlix) >u= 2m < j^ix) >u + < A(x) >u ■ (9) 

The lattice definitions of the various (gauge invariant) operators in (9) are only 
unique up to terms vanishing in the naive continuum limit. Hence the choice of the 
left derivative in (9) is only a matter of convenience. It appears naturally, e.g., when 
Wilson fermions are considered. 

Let 0{x) stand for any of the operators appearing in (9). Then < 0{x) >u has 
the following formal expansion in the gauge potentials 

< 0[X) >^= E Tf;^X"^^\^^^^^^---^^n)Al\{x^)---A^:{Xn) , (10) 

n>2 ■ {xi},{tii},{ai} 
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where, because of the assumed symmetries of (1), the sum over n starts with n = 

2. The correlation functions r^^.^."/l„""" ( symmetric under the 
exchange of any pair of collective labels {xi,fii,ai). Defining the Fourier transform 

of t''^^^^^'„'°'"{x\xi,X2, ■■■,Xn)hy 

2L (27r)4 11 (27r)4 ^m-Mn Wl^i' ' '^^"J ' 

(11) 

where, by translational invariance, 

n 

f . . ., = 5(g - ^ h)rfj:';,-r-{h, ■ ■ -, K) , (12) 

i=l 

the Ward identity (9) translates as follows to momentum space, 

-^9;.%;^%„(^ir--,M = 2mf^--:;,J/.i,---,fc„) 

+ f(t.).X"(^i'---'M , (13) 

where 

= 2 - sm . (14) 

As we shall see further below, gauge invariance implies that every term in eq. (13) 
possesses a finite continuum limit. Any anomalous contribution to the divergence 
of the axial vector current must therefore have its origin in the naively irrelevant A 
term. If its continuum limit is universal, then so is the anomaly. 

We now state the central theorem of this paper and then give the proof: 

3. Theorem 

Any lattice discretization of the action (1) with the following properties: a) S 
has the correct continuum limit; b) S is gauge invariant; c) The Dirac operator is 
local and d) The free Dirac operator is invertible for nonvanishing momenta (no 
doublers), reproduces the axial anomaly in the continuum limit. 

Let us consider first in turn the consequences of b) to d). 

i) Gauge Invariance 

Gauge invariance tells us that if (9(74,-0,7/;) is a gauge invariant operator, then its 
external field expectation value satisfies 

< 0(A", 0, ^) >A"=< 0{A, ^, ^) >A , (15) 

where is the gauge transformed potential. On the lattice the variation of the 
gauge potentials induced by an infinitessimal gauge transformation is given by P 

^Al^x) = [gfatcA'Jx) - M-\gaA,{x))d;y{x) , (16) 



where fabc are the structure constants of SU{3), is the dimensioned right lattice 
derivative, and the matrix M is given by 



1 _ „-i<j){x) 

with (p a colour matrix in the adjoint representation. A summation over repeated 
indices is understood. Up to 0{(f)'^) we have that 

M-^(0) = l + ^0-^(0^) + W^) . (18) 

Performing the transformation (16) in (10), one is led to an infinite set of gauge 
identities. Only the following two will however be required for the proof of the 
theorem: 

^^r^Ji^^'^HxIxix^) = , (19) 



^IJ-l^iil}i2lJ'3 ^(^1^12^23^3) 9^X1X2 faia2C^li2lii ^ {•^{•^'iXs) 

- l^agUa2cd;:jx,.,TfJ;;^^ix\x2x,) 

+ (2^3), (20) 

where the left-derivatives acts on Xi. Written in momentum space (19) takes the 
form 

CK),.i'^^'.l7''{ki,h) = o, (21) 

where * denotes complex conjugation, and where is defined in an analogous way to 
(14). Because of Bose symmetry a corresponding statement holds for (ki)^^ replaced 
by (A;2)^2- The second identity (20) reads as follows in momentum space, 

^Ckl),Jf^l7;rikl, k2, ks) = 9fa.a2cr'SlTiki + h, ks) 

- ^ay,.a2c{kl),2^T,l^h + k,M) 

+ (2^3). (22) 



ii) Locality of the Dime operator 

The coefficient functions D^"^"Jj'„'""(x, y\xi ■ ■ ■ x„) in (2) are assumed to vanish at 
least exponentially fast for large separations between any pair of lattice sites with a 
decay constant of the order of the inverse lattice spacing. This is trivially the case for 
Wilson fermions, where the Dirac operator connects only neighbouring lattice sites. 
It also holds for Ginsparg- Wilson fermions, for sufficiently smooth background gauge 
field configurations . As a consequence the correlation functions f are 
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analytic functions of the momenta around vanishing momenta. This fact, combined 
with eq. (21) (following from gauge invariance), implies that 

rif(SrH^i'^2) = 0, (23) 

where Ti denotes the Taylor expansion around zero momenta up to first order. 

Consider next eq.(22) evaluated for k2 = = 0, and small ki. Because of (23) 
the RHS is of 0{kl). It therefore follows that 

fi';'i:rr='(o,o,o) = o. (24) 

Relations (23) and (24) are weaker than those obtained from gauge invariance in 
the case of QED, where an analogous statement to (23) holds for vertex functions 
involving an arbitrary number of external photon lines. They nevertheless suffice to 
prove the theorem. 

iii) The free propagator D^'^\p)^^ is free of doublers 

As a consequence the lattice power counting rules and the Reisz theorem ^ 
applies. This theorem allows us to take the naive continuum limit of a Feynman 
integral, if the lattice degrees of divergence of all Zimmermann spaces are negative. 

4. Proof of the theorem 

We are now ready to prove the theorem. Our attention will be focused on the 
"irrelevant" contribution of the A-term in (13), which should generate the anomaly 
in the continuum limit. 

Consider the axial vector Ward identity (13). By power counting the ultraviolet 
lattice degree of divergence (LDD) of T^^ and Fs is given by 3 — t?,, where n is the 
number of external gluon fields, while the LDD of F^^-* is 4 — n. Hence the LDD 
of Feynman integrals contributing to F^^^ is negative for graphs involving more 
than 4 external gauge fields. Their contributions thus vanish by the Reisz theorem 
PI in the continuum limit, since A is an irrelevant operator. We therefore only 
need to consider the correlation functions T^f^l^^'^'^iki, k2), f AnM2M3^"^(^i) ^2, ^3), and 
f |itM2M3M?°*(^ii ^2, ^3, ^4) with LDD = 2,1 and 0, respectively. Let us decompose 
these vertex functions as follows 

fitir(fci'^2) = {l-T2)T^,'^lrHki,k2) 

+ T.tl'^lT'iki,^) + {T,-T,)V^^^ll-'^-{k,,k2) , (25) 

i'iVz7^'ik,,k2, k,) = (1 - Tof(tirr'(A;i, k,) + Tof(ti:j."rHfci, ^3) 

+ (Ti-To)f(t)^--^(A:i,fc2,A:3) , (26) 

and 

riti:Mrr'(^i' ^2, k,) = (i - To)f(t)^;,x'^^(fci, k^, k,, k,) 
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+ Torl%Z7:r(ki,k2,h,h) , (27) 

where T„ denotes the Taylor expansion around vanishing momenta up to n'th or- 
der. By the Reisz theorem the first term appearing on the RHS of (25) to (27), 
respectively, vanishes in the continuum limit, since it has negative LDD and A is an 
irrelevant operator. The respective second terms in (25) and (26) vanish by gauge 
invariance (cf. eqs (23) and (24)). Hence we conclude that 

}}^/i'^lTHh,h)^^miT,-T^)^^^^^^^^ , (28) 

li^n^™"'(^i' ^2, h) = limfTi - ro)f (^--'^^(^^^ j,^^ ^ (29) 

li^o^M™r/"(^i' = li^o^o^™"^^^^ • (3°) 

As we now show the limits on the RHS are universal, i.e., they do not depend on 
the specific form of the A operator in (9). 

Consider first the RHS of (25). Applying the operation {T2—T1) to the axial Ward 
identity (13) with n — 2, and making use of (23), following from gauge invariance, 
we conclude that 

(T, - Ti)f(,'^:^"^(fci, h) = -2m{T, - Ti)f-- ,(/.!, h) . (31) 

The expression appearing on the RHS has negative LDD; hence its continuum limit 
is given, according to the Reisz theorem, by applying the (T2 — Ti) operation to 
the integrand of the corresponding continuum Feynman integral, i.e. the triangle 
diagram. One then readily finds that (28) is given by 

lim f (fi^^'^^/ci, k^) = -£^5a,a,e,,,,.,{hUh), . (32) 

Consider next the RHS of (29). Applying the (Ti— Tq) operation to the Ward identity 
(13) with n — 3, and making use of (24) (following again from gauge invariance) 
one obtains 

(Ti - To)f(t)^J,^^«3(^^^ f,^^ ^.3) ^ _2m{T, - To)^}^^^, h, h) . (33) 

Since the LDD of the RHS is negative, its continuum limit is given by applying 
the (Ti — To) operation to the integrand of the continuum box Feynman diagram, 
involving three vector vertices and a 75 insertion. After some work one finds that 

where q = ki + k2 + k^. Finally, consider (30). Applying Tq to the Ward identity 
(13) we obtain 

li«?,^of (ti~-(fci, k2, ks, k,) = -2mlim rof~3V,(A;i, k,, k^, k,) . (35) 
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Since the LDD of the vertex function on the RHS is again negative, the hmit a — > 
can be calculated by applying Tq to the continuum pentagon graph involving a 75 
insertion. A simple calculation shows that (35) is given by 

^^/lZ',7:riki,h,h,h) = |^[eM,;.,,3M4rr(T»^T-r-T''^) +perm.] (36) 

where T" are the SU{3) generators in the fundamental representation. The (con- 
tinuum) correlation functions in coordinate space, corresponding to (32), (34) and 
(36), are obtained by performing the Fourier transformation (11), where the lim- 
its of integration now extend to infinity, and with t^i^\"'}'f^'^'^ defined in (12). The 
anomalous contribution is then given by (10) with O — > A, where the sum over the 
coordinates are replaced by integrals. Only the two and three gluon vertex functions 
actually contribute to the anomaly, as can be shown by making use of the Jacobi 
identity. One then easily verifies that the anomalous contribution to the (euclidean) 
axial vector Ward identity takes the well known form 

^^""^ ^ - J^e^-ApF;.(x)F,;(x) (37) 

where 

F;,(x) = d,Al{x) - d^Alix) - gUcK{x)At{x) (38) 

is the non-abelian field strength tensor. 

Finally, we point out that not only f ^'^.^."/l^"" possesses a finite limit, but also 
the other two vertex functions appearing in (13). The reason is, that the latter have 
LDD = 1 and for n = 2 and n — 3, respectively, and negative LDD for n > 4. 
Gauge invariance, as embodied in the statements (23) and (24), allows us to replace 
the n = 2 and n = 3 vertex functions by the corresponding Taylor subtracted forms 
with negative LDD. These possess a finite continuum limit. 

Concluding, we have shown, that any choice of gauge invariant lattice discretiza- 
tion of the action of QCD which is local in a more general sense and free of doublers 
will reproduce the correct axial anomaly. This anomaly, which arises from a naively 
irrelevant term in the Ward identity, is determined from the appropriate Taylor 
terms in the expansion of the triangle, box and pentagon continuum Feynman dia- 
grams around vanishing momenta. Although we have explicitely referred to QCD, 
it is evident that the proof holds for any SU{N) gauge theory. 
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